Classical-interference analog of quantum fluctuations for bound-state soliton pairs 
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Quantum photon-number fluctuation and correlation of bound soliton pairs in mode-locked fiber 
lasers are studied based on the complex Ginzburg-Landau equation model. We find that, depending 
on their phase difference, the total photon-number noise of the bound soliton pair can be larger or 
smaller than that of a single soliton and the two solitons in the soliton pairs are with positive or 
negative photon-number correlation, correspondingly. It is predicted for the first time that out-of- 
phase soliton pairs can exhibit less noises due to negative correlation. 
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Quantum solitons have attracted a great deal of re- 
search interest in the contexts of nonlinear quantum op- 
tics, condensed-matter physics, and quantum informa- 
tion science due to their remarkable nonclassical prop- 
erties. In particular, quantum solitons in optical fibers 
largely resemble their classical counterparts, but with ad- 
ditional quantum fluctuations around the mean fields. It 
has been possible to achieve squeezing through quantum 
solitons in optical fibers, 0, 0, 0, an d they may also 
serve as a new platform for quantum information appli- 
cations HH0- 

Quantum solitons are macroscopic optical wave pack- 
ets which offer a tcstbed for quantum optics and quantum 
field theories. For the quantum nonlinear Schrodinger 
equation (NLSE), exact soliton states can be constructed 
as combinations of eigenstates of the Hamiltonian of the 
one-dimensional Bose gas with 5-like (contact) interac- 
tion through the Bethe ansatz method In the large 
photon number limit, which corresponds to the usual op- 
tical solitons generated by lasers, the many photon wave 
function of the quantum soliton is well approximated 
by a single-photon wave function (the Hartree approx- 
imation) [j| . Linearization around such a soliton j^, [j| 
successfully explains experimental observations of quan- 
tum fluctuations for temporal fiber solitons, provided 
that optical loss and higher-order effects are negligible 
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It is well known that the force between adjacent soli- 
tons in the NLSE model is attractive or repulsive, de- 
pending on the phase difference between them |I5|. Sta- 
tionary bound soliton states in this conservative model 
do not exist. Formation of effectively stable double-, 
triple-, and multi-soliton bound states was predicted in 
models based on the complex Ginzburg-Landau equation 
(CGLE) 0, 0, , and observed experimentally in var- 
ious passively mode-locked fiber lasers 0, |2jJ, E| • The 
separation between the solitons in these bound states 
are "quantized", taking a set of discrete values. The 
amplitude noise in triplet bound states generated by a 
stretched-pulse ytterbium-doped double-clad fiber laser 
was observed to be reduced compared to the single soli- 
ton pulse |23. It is an issue of straightforward interest 
to study the noise of these bound solitons, and to un- 



derstand why the mode-locked fiber lasers operate more 
stably in the bound-state regime. 

The passively mode-locked fiber lasers are quite accu- 
rately described by the cubic-quintic CGLE. In a normal- 
ized form, the equation is 



iU z + (D/2)U tt + \U\ 2 U 



= iSU + ie\U\ 2 U + if3U tt 
+ i^\U\ 4 U -is\U\ 4 U, (1) 



where U is the local amplitude of the electromagnetic 
wave, z is the propagation distance, t is the retarded 
time, and D = +1 and —I correspond, respectively, to 
the anomalous and normal dispersion. Besides the group- 
velocity dispersion (GVD) and the Kerr effect, which are 
accounted for by conservative terms on the left-hand side 
of Eq. (TIJ, the model also includes the quintic correction 
to the Kerr nonlinearity, through the coefficient u, and 
non-conservative terms. The coefficients S, e, fi, and (3 
account for the linear, cubic, and quintic loss or gain, and 
spectral filtering, respectively. 

In the CGLE model, with suitable parameters degener- 
ate bound-state soliton pairs are known to exist through 
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FIG. 1: (A) The photon-number correlation parameter, C12, 
for the bound-state soliton pairs with different relative phases. 
(B) The contour plot of the classical solution for the bound 
state. All the results are presented for D = 1, 8 = —0.01, 
e = 1.8, (3 = 0.5, ft = -0.05, and v = in the CGLE Q. 
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the balance between the gain and loss, in the form 
ilH, U(z,t) =U (z,t + p)e- ie / 2 + U (z,t-p)e ie / 2 , 
where Uq is a single soliton solution, and p and 9 are 
the separation and phase difference between the soli- 
tons. In this Letter, we focus on the consideration of 
three fundamentally different cases, corresponding to the 
bound states with the same separation and amplitude, 
and 6 = 0, tt/2, and tt (the in-phase, orthogonal, and 
out-of-phase pair), respectively. 

We compute the quantum fluctuations of these soli- 
ton pairs by dint of a numerically implemented back- 
propagation method |2.'"'J| . which may be summarized as 
follows. First of all, we replace the classical function 
U(z,t) in Eq. Q by the quantum-field operator vari- 
able, U(z, t), which satisfies the equal-coordinate Bosonic 
commutation relations. Next, the equation is linearized 
around the classical solution through the substitution of 
U(z,t) = Uo(z,t) + u(z,t), assuming large photon num- 
bers in the solitons. Then, a zero-mean additional noise 
operator, n(z,t), is introduced to make the quantum 
perturbation fields in the linearized equation satisfy the 
Bosonic communication relations (see Ref. y for more de- 
tails). By imposing suitable correlation functions for the 
noise operator, the minimum quantum noise in the con- 
sidered nonconservative model is introduced. Therefore 
the results presented here represent a lower limit required 
by the fundamental principles of quantum mechanics. 

Figure ^ shows the photon- number correlation param- 
eter for the two solitons in the bound soliton pair, which 
is defined as 



(A) 8 = 



c 



12 



ANiAN 2 :) 
(AiVf )(AiVf) ' 



Here, the colons stand for the normal ordering of the 
operators and ANi^ are perturbations of the photon- 
number operators for the two solitons, which are num- 
bered (1,2) according to their position in the time do- 
main. Initially, the two solitons are assumed to be un- 
correlated, with fluctuations around each soliton obeying 
the coherent-state statistics. For the in-phase pair, the 
photon-number correlation between the solitons gradu- 
ally increases to positive values and eventually saturates 
around C\i = 0.36. But for the out-of-phase pair, C12 
gradually decreases to negative values and then saturates 
too. In between, the correlation parameter for the case 
of 6 = tt/2 remains close to zero as long as the computa- 
tion is run. For the former two cases, the saturation of 
the photon-number correlation parameter is due to the 
nonconservative effects in the CGLE model. 

To further demonstrate the behavior difference of the 
photon-number correlation for soliton pairs with differ- 
ent relative phases, in Fig. Elwe display the time-domain 
photon-number correlation patterns for them. The plot- 
ted correlation coefficients, , are defined through the 




(C) 8 = 11 



FIG. 2: The time-domain photon-number correlation pat- 
terns, rjij, for the bound soliton pairs with different relative 
phases, after the normalized propagation distance z = 0.4. 
(A): 8 = 0, (B): 8 = tt/2, and (C): = it. The time-division 
length At = 0.3 



normally-ordered covariance, 
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where Ahj is the photon-number fluctuation in the j-th 
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FIG. 3: Comparison of the total photon- number fluctuations 
in the bound soliton pairs with different relative phases (three 
solid lines), and around the single soliton (the dash line). 

time slot Atj, 

Ahj = / dt[U (z,t)u [ (z,t) + Uo(z,t)u(z,t)]. 

J Atj 

Here the integral is taken over the given time slot, with 
the same time-division length At. Clearly, in Fig. [21 
(A) one can see that there is a strong posiiwe-correlation 
band connecting the quantum correlation patterns of the 
bound solitons when they are in phase, 9 = 0. In Fig. 
121(C) there exists a negaiive-correlation pattern between 
two solitons for the out-of-phase case, 9 = tt. Moreover, 
for the case of 9 = tt/2, in Fig. 121(B), the correlation pat- 



terns of bound solitons are almost isolated. In classical 
physics, in-phase and out-of-phase fields will lead respec- 
tively to the constructive and destructive interference. 
Here we observe a similar effect for the quantum noises. 
What is more important, in Fig. |31we compute the total 
photon number noise of the bound soliton pair and com- 
pare it to the case of a single soliton (these results are 
amenable to straightforward experimental verification). 
As one may expect, the photon-number noise of the in- 
phase pair is larger than that for the single soliton, which 
may be explained as the fluctuation enhancement due to 
constructive interference. On the other hand, the noise 
is reduced for the case of out-of-phase pair as the result 
of destructive interference. The orthogonal soliton pair 
with 9 = tt/2 may be viewed, in the first approxima- 
tion, as independent two single solitons, which explains 
why it features almost the same noise level as the single 
soliton, even though small oscillation of the noise level 
originated from the residual interaction between the two 
solitons can still be seen. 

In conclusion, we have presented theoretical results on 
the photon-number correlation and total photon-number 
noise for bound-state soliton pairs in the model of com- 
plex cubic-quintic Ginzburg-Landau equation. The cases 
of the in-phase, orthogonal, and out-of-phase soliton 
pairs have been considered in detail. We conclude that 
the interference of the quantum fluctuations in the soli- 
ton pair is constructive or destructive depending on the 
classical relative phase of the solitons. An important 
consequence of the results is that the operation regime of 
the fiber laser should be more stable when it is based on 
the out-of-phase soliton pairs. 
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